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Summary. This is an introduction to the study of the in-medium behavior of 
quarkonia and its application to the quark-gluon plasma search in high energy nu- 
clear collisions. 



1.1 What are quarkonia? 

The bound states of a heavy quark and its anti-quark which are stable with respect 
to strong decay into open charm or bottom are collectively called quarkonia. We 
denote by Q either of the heavy quarks, charm (c) or bottom (b); the corresponding 
bound states are known as charmonia or bottomonia, respectively. 

Among the vector (spin-one) charmonium states, the lightest (ground state) is 
the famous J/tp; the excited states are the Xc and the tp' . For the bottom quark, the 
lightest quarkonium is the T, while the excited states include the xt, ^' , x'b s-nd the 
T". The stability of the cc/bb quarkonium states implies that their masses satisfy 
Mcc < 2Md and A/^j, < 2Mb, where D — cu and B — bu are the corresponding 
"open" mesons. A specific characteristic of quarkonia is their small size. While the 
typical hadron radius is ~ 1 fm, the radii of charmonia and bottomonia range from 
0.1 — 0.3 fm, as we shall see. 

Since c and b quarks are very heavy {rriQ <C Aqcd ~ 200 MeV), the binding of 
the QQ system may be treated non-relativistically. The governing equation is the 
non-relativistic Schrodinger equation, 

- — {V'(r) + V{r)}<I',(r) = (M, - 2m)'j^,{r), (1.1) 
m 

where \['{r) denotes the wavefunction of the system, r the quark-antiquark sepa- 
ration, and m the quark mas tQ. Since Eq. IlTJ 

is a nonrelativistic description of 
the binding, the total rest mass must be subtracted from the masses Mi of the 



^ We work in the center-of-mass system, with a reduced mass m/2, so that we have 
— V^/m instead of the usual — V^/2m. 
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Xb 
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mass (GeV) 


3.10 


3.53 


3.68 


9.46 


9.99 


10.02 


10.36 


10.36 


AE (GeV) 


0.64 


0.20 


0.05 


1.10 


0.67 


0.54 


0.31 


0.20 


radius (fm) 


0.25 


0.36 


0.45 


0.14 


0.22 


0.28 


0.34 


0.39 



Table 1.1. Masses, binding energies and radii of lowest cc and bb bound states [2]. 



bound states. Once we find the eigenvalues Mi of the system, we can also define the 
"binding energy" AE of each quarkonium state, AE = 2Md,b — Mi. 

Lattice and spectroscopic studies suggest for the potential V{x) the form [1] 

VM^ar--, (1.2) 
r 

generally known as the "Cornell potential". It is spherically symmetric, and consists 
of two parts. The linearly rising part represents the confining force, given in terms of 
the string tension a; lattice studies put its value at around 0.2 (GeV)^. The second 
part is an effective Coulomb potential, including transverse string oscillations; string 
theory suggests a = 7r/12. 

Having solved the Schrodinger equation, we may determine the bound-state radii 
through 

^ Jd3r\^,{r)\^ ■ ^ ' 

A fair estimate can already be obtained by means of a semi-classical formulation. 
The energy of the system is then given by 

E=^ + V{r), (1.4) 
m 

and from the uncertainty relation we have pr ~ c; the constant c can be fixed by 
requiring the correct J/tp mass, giving c ~ 1.25. Minimizing the energy determines 
the lowest bound state radius ro, 

cr-f- — = 3. 1.5 

With a ~ 0.2(GeV)^ and a ~ 7r/12, together with rric — 1.3 GeV, we obtain a 
J/ip size {QQ separation, i.e., twice the radius) of about 0.5 fm. For the a' = 
value, we have ro ~ {l/ma)^^^ ~ 0.3 fm; on the other hand, for a = 0, we get 
ro ~ (1/ma) ~ 0.6 fm. We thus see that a major contribution to the radius comes 
from the string tension. At T = 0, the radius of the J/tp is thus to a considerable 
extent still determined by the confining part of the potential. We summarize some 
of the characteristics of the spin-averaged quarkonia in Table [01 [2]. 



Next, we turn to the question of the dissociation and decay of quarkonia. We 
have already noted that these mesons cannot decay via strong channels because 
their masses are smaller than the open thresholds. It is also known that quarkonia 
do not dissociate significantly in nuclear collisions; we shall discuss this in greater 
detail in section [L3] 
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How then do quarkonia dissociate? Three mechanisms have been identified, cor- 
responding to the behavior for T — 0, < T < Tc and T > Tc, where Tc is the 
critical temperature of deconfinement. We shall consider each of them in turn. 



1.1.1 String-breaking 

The potential in Eg. I1.2l is correct only in the limit Mq ^ cx). If light quarks exist 
in the theory, then the string connecting the heavy quarks can break as soon as the 
overall energy in the system is greater than 2Mn or 2Mb, depending on whether 
Q = c or b (Fig. . Light quark-antiquark pairs appear at the broken ends of the 
string, and new "heavy-light" Qq or qQ mesons are formed. This behavior of the 
quark potential has been observed in lattice studies with dynamical quarks, as we 
shall show below. 



Q Q 



Q Q 



4 ► 



49* 



Fig. 1.1. Cartoon of string-breaking 



We may estimate the string-breaking energy Fq. For the charm quark, Fo = 
2{Md - rric) ~ 1.2 GeV, while for the bottom quark, Fo = 2{Mb - mt) ~ 1.2 GeV. 
From this, we deduce ro = (1.2 GeV) /a ~ 1.5 fm. That this value is the same for 
both quark species leads us to conclude that the energy required for string breaking 
is a property of the vacuum itself, as a medium at T = 0, containing virtual qq 
pairs which are brought on-shell by the field between the heavy quarks. The effect 
of string breaking on the Cornell potential is shown in Fig. 11.21 



1.1.2 Recombination 

In nuclear collisions not sufficiently energetic to create a quark-gluon plasma, there 
will nevertheless be abundant hadron production. These newly formed light hadrons 
can through a switch in bonding (recombination) turn a QQ meson into two heavy- 
light mesons. This mechanism is schematically depicted in Fig. 11.31 when two or 
more hadrons overlap, their quarks can recouple form new pairs. 

If the temperature is increased, the hadron density also increases, and this in 
turn increases the recombination probability. As a consequence, the distance up to 
which the heavy quarks still bind also becomes shorter, and the potential will break 
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Fig. 1.2. String breaking energy as function of QQ separation 




Fig. 1.3. A schematic view of recombination 
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Fig. 1.4. Schematic dependence of the string-breaking radius with temperature 



earher (see Fig. 11.41 We thus have something hke "effective screening" , even though 
all color charges are still bound. 

What happens as we get close to Tc? The density of produced hadrons will then 
increase strongly, and lattice studies show that in accord with our picture, both the 
free energy and the string-breaking radius tt decrease rapidly near Tc, as shown in 
Fig.O 

1.1.3 Color screening 

Above T — Tc, we have a medium of unbound color charges, and an entirely different 
mechanism takes over. At all temperatures T above zero, quarks and gluons are 
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Fig. 1.5. Lattice results for free energy and screening radius as function of T. 
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Fig. 1.6. Quarkonium dissociation as "thermometer" for the quark-gluon plasma. 



screened, just as electric charges experience Debye screening in an electromagnetic 
plasma. This screening occurs with a characteristic radius, which we denote by 
td. It decreases with increasing temperature, as the medium increases in density. 
Deconfinement is expected to occur when this radius becomes comparable to the 
average hadron size of 1 fm. Then a given quark can no longer see its former partner 
in a hadron; instead, it sees many other quarks and antiquarks and therefore can 
move around freely, without encountering any confinement limit, since it is never 1 
fm away from an antiquark. 

We would like to use the behavior of the J/tp to probe if a quark-gluon plasma 
was formed in the collision [5]. However, as we have seen, the J ftp and its heavier 
counterparts have smaller radii than the usual hadrons. Thus, charmonia and bot- 
tomonia may be expected to survive beyond the QGP phase transition up to some 
higher temperature, at which they will become dissociated. Thus, if we know their 
sizes as well as the behaviour of rn as a function of T, we can use their dissociation 
points to determine the temperature and the energy density e of the QGP medium 
[i], as illustrated in Fig. 11.61 
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1.2 Studying charmonium dissociation 

We now turn to the question of how to determine quantitatively the quarkonium 
dissociation points in a quark-gluon plasma. Two different approaches were used to 
address this problem. 

• Solve the Schrodinger equation with a temperature-dependent potential V{r, T), 
or 

• calculate the quarkonium spectrum directly in finite temperature lattice QCD. 
We shall look at each of these approaches in turn. 

1.2.1 Potential models for quarkonium dissociation 
The Schwinger model [Sj 

One generalizes the Cornell potential, eq. (|1.2|l , to non-zero temperature in the form 

V{r,T) = arl^^^^-^\ ~-e~^"'. (1.6) 
I fir ) r 

The T-dependence of the above equation is in the "screening mass" fJ.{T) = l/roiT). 
Eq. (|1.6|l gives the correct zero-temperature limit, eq. (|1.2|) . for jJ.{T) as T 0. 

To determine the dissociation points, one solves the Schrodinger equation and 
determines the bound-state energies Mi{fi). With increasing temperature, the bound 
state i disappears at some fi = fn. One then uses the temperature dependence of 
the screening mass from lattice estimates, fJ,{T) ~ 4 T, to determine the Ti. The 
result of this model is 

• the ip' and Xc become dissociated around T ~Tc, 

• the J/ip survives up to about T ~ 1.2Tc. 

In both cases, at the dissociation point the binding energy vanishes, while the 
binding radius diverges. 

Lattice potential models O E H [9l [lO] 

Alternatively, one may use lattice results for the temperature dependence of the 
potential felt by a static quark-antiquark pair to determine the needed potential. 
The static QQ studies start from the partition function Z, which is related to the 
free energy hy Z — exp(— /3_F); this in turn gives the thermodynamic potentials 

F = U~TS 

U{r,T) = F{r,T)-T(^^^^^^ (1.7) 

Assuming that the internal energy U (r, T) provides the temperature dependence of 
the heavy quark potential, we we use results from Nf = 2 lattice QCD and solve 
the Schrodinger equation. The results obtained from such studies indicate that 
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Fig. 1.7. Schematic view of lattice results for charmonium dissociation. 

• the and Xc are dissociated around a temperature T ~ l.lTc, 

• J/ip survives up to a temperature T ~ 2 Tc. 

Comparing these results to the ones from the Schwinger model, we see that while 
there is agreement in the case of the higher excited states, lattice potential models 
predict a considerably higher dissociation temperature for the J/'ip. The reason for 
this is that the internal energy U{r, T) leads to much stronger binding than the 
Schwinger model potential. 

It should be noted here that there still is some ambiguity as to whether U or F 
is the correct potential to be used in the Schrodinger equation. Hence there exist 
approaches with potentials of the form all +{l — a)F, with < a < 1. Such potentials 
tend to reduce binding and lower the dissociation temperature as a is decreased. 

1.2.2 Lattice studies of charmonium survival 

The ideal way to resolve the above ambiguity would be to calculate the cc spectrum 
directly on the lattice, and this is indeed what lattice studies aim to do [11II12|[T51|14I 
1151116) . More specifically, they calculate the cc spectrum a{u!,T) in the appropriate 
quantum channel, as a function of the temperature T and the cc energy uj. Bound- 
states show up as resonances in a plot of a versus ui. By performing simulations 
at different temperatures, one can determine the temperature at which a particular 
peak disappears i.e. a bound-state dissolves. A schematic illustration is shown in 
Fig. 11.71 The results presently indicate that 

• Xc is dissociated for T > l.lTc. 

• J/ip persists upto 1.5 < T/Tc < 2.3. 

Thus, on the basis of lattice studies, the following picture emerges: The higher 
excited states dissociate around T — Tc, while the J/^ survives up to much higher 
temperature, in accord with the potential model studies based on the internal energy 



There is, however, a caveat to these calculations. The discretization introduced 
by the lattice limits the resolution of the peak. Lattice methods are thus useful 
in determining the position and to some extent the amplitude of the peaks, but 
determining the peak widths remains a challenge, nor is it easy to study the spectrum 
in the continuum region (tj > 4 GeV). 



U{r,T). 
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1.3 Dynamics of quarkonium dissociation 

We have seen in the previous discussion that the J/V', the vector ground state 
of charmonium family, is very tightly bound. Its binding energy i.e. the energy 
difference between J/ijj mass and open charm threshold, AEj/^, is considerably 
larger than the typical non-perturbative hadronic scale Aqcd, 

AEj/^ = 2Md - Mj/^ = 0.6 GeV > Aqcd. ~ 0.2 GeV (1.8) 

Consequently the size of J/ip is much smaller than that of typical hadron, 

rj/^ ~ 0.25 fm < Aq^.^ ~ 1 fm (1.9) 

We now want to consider by what kind of dynamical interaction such a state can be 
dissociated. Because of the small spatial size, the J/ip can only be resolved by a suffi- 
ciently hard probe. Moreover, because of its high binding energy, only a sufficiently 
energetic projectile can break the binding. The previous study of global medium 
effects had led to the conclusion that only a hot deconfined medium, consisting of 
colored quarks and gluons, is capable of dissociating the charmonium vector ground 
state. We now want to study this on a microscopic level. 

In a deconfined medium, the constituents are unbound partons, whereas in a 
confined medium the constituents are hadrons. Such thermal hadrons are incapable 
of causing coUisional dissociation of J/tp. Let us illustrate this point. 

Consider the collision of a J/tp with a normal hadrorn. Because of the small 
characteristic J/tp size, only a hard partonic constituent of the hadron can see the 
J/tp and interact with it. In other words, J/tp collisions with ordinary hadrons probe 
the local partonic structure of these 'light' hadrons, not their global aspects such as 
mass, size, or overall momentum. The parton nature of the interaction is illustrated 
in Fig. [TSI 



JA|f 




Fig. 1.8. Schematic view of the interaction of a normal hadron with a J/V"- 



To see the effect of this more quantitatively |17) . we take an ideal pion gas as 
the confined medium. The momentum distribution of the pious at a temperature T 
follows f{p) ~ exp(— |p|/r), giving the pious an average momentum (p) ~ ST. Now 
the gluon momentum distribution inside a hadron as determined by deep inelastic 
lepton- hadron scattering is given by the parton distribution function g{x); here 
X — 2kg /y'^, with kg for the gluon momentum, so that x may be thought of as the 
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fraction of the incident liadron's momentum carried by tlie gluon. For tiie pionic 
gluon it takes tlie form 

g{x) ~ {l-xf. (1.10) 
The resulting average gluon momentum in the hadron thus becomes 

\ lodx X g(x) 

{kg)=Ph--^ — . (1.11) 

Jo dx g{x) 

With eq. (|1.10|l and ph ~ 3T for the momentum of the incident hadron, we obtain 

{fc«). = ^ = ^<0.1GeV, (1.12) 

where we have assumed T < 175 MeV for the temperature of the hadronic medium. 
Thus gluons bound inside the hadronic constituents of confined matter are much 
too soft to cause the dissociation of a J/tp. 

On the other hand, in a deconfined medium, such as an ideal QGP, the gluons 
are free and distributed according to a thermal distribution f{kg) ~ exp(— fcj/T), 
which gives 

(k,) ~ 3T (1.13) 

so that for T > 1.2 Tc — 0.63 GeV, the gluons are hard enough to overcome the 
J/^ binding. 

We have thus noticed that deconfinement results in a hardening of the relevant 
gluon momentum distributions. More generally speaking, the onset of deconfinement 
will lead to parton distribution functions which are different from those for free 
hadrons, as determined by DIS experiments. Since hard gluons are needed to resolve 
and dissociate J/^p's , one can use J/^ps to probe the in-medium gluon hardness and 
hence the confinement status of the medium. 

This qualitative picture can be made quantitative by short distance QCD calcu- 
lations [171 [18] . One has to calculate first the cross section for gluon dissociation of 
J/V', a QCD analogue of the photo-effect. This can be carried out using the operator 
product expansion, which is essentially a multipole expansion for the charmonium 
quark-antiquark system. Fig. 11.81 shows the relevant diagram for the calculation of 
inelastic J/ip-haAron cross section. The upper part of the figure corresponds to J/ip 
dissociation by gluon interaction. The cross section for this process has the form 

1 jk/AE^ - 1)«/^ 

with AEj/^ = 2Md — Mj/^. The corresponding cross section for the hadron dis- 
sociation is obtained by convoluting this gluon dissociation cross section with the 
gluon distribution function g{x) of the incident hadron. For J/7/;-meson interactions, 
this leads to the form ^ 

i^h-j/ii — o-goom(l r-)^'^ (1-15) 



with A ~ (s - Ml)/M^ and Aq ~ {Mn ^ AE^), where ^/s is the CMS energy of 
the J/t/j- hadron system. Here cTgoom — Tj^^ ~ 2 mb is the geometric J/ip cross 
section and Mh denotes the mass of the incident meson. Fig. 11.91 compares the two 
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Fig. 1.9. Cross sections for J/ip dissociation by gluons vs. pions 

dissociation cross sections, J/tp dissociation by gluons (gluo-efFect) and by pions, as a 
function of projectile momentum k incident on stationary J/V", as given by eqs. p.l4[> 
and (|1.15|) . The gluon cross section shows the typical photo-effect form, vanishing 
until the gluon momentum kg reaches the binding energy AEj/^; it peaks just a 
little later {Xg ~ tj/^) and then vanishes again when sufficiently hard gluons just 
pass through the (comparatively larger) charmonium bound states (Ag ^ tj/^). In 
contrast, the J/i/)-hadron inelastic cross section remains negligibly small until rather 
high hadron momenta (3-4 GeV). In a thermal medium such momenta corresponds 
to temperatures of more than 1 GeV. In other words, in a confined medium in the 
temperature range of the order of a few hundred MeV the J /tj) should survive, but 
it should become dissociated in a hot deconfined medium. Confined media in the 
temperature range of few hundred MeV are thus essentially transparent to a J/i), 
while a deconfined medium of the same temperature is opaque to J ftp's and very 
efflciently dissociates them. 



1.4 Quarkonium production in nuclear collisions 

The aim of ultra-relativistic nuclear collisions is to study color deconfinement and the 
resulting quark-gluon plasma in the laboratory. We want to use quarkonia produced 
in the collision as a probe to study the medium produced in the collision. Both 
the quarkonium states and the medium to be probed require a 'finite formation 
time', so we have to look at the evolution aspects in both cases. Let us first consider 
the issue of charmonium production in hadron-hadron collisions and then turn to 
nuclear targets. 

1.4.1 Quarkonium production in hadronic collisions 

Quarkonium production in hadron-hadron collisions occurs in three stages. The first 
stage is the production of cc pair. Because of the large quark mass {rric ~ 1.3 GeV ) 
this process can be treated as a hard process and is well described by perturbative 
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QCD. A parton from the projectile interacts with one from the target; the (non- 
perturbative) parton distributions within the hadrons are determined empirically in 
other reactions, e.g., by deep inelastic lepton-hadron scattering. At high energy the 
process of cc production dominantly occurs by gluon fusion, gg cc (see Fig. ll.lOjl . 



Fig. 1.10. Lowest-order Feynman diagram for cc production through gluon fusion. 



The cc in general is in a color octet state. It has to neutralize its color in order 
to leave the interaction zone and form a physical resonance like J/ip or ip' . In the 
second stage, color neutralization occurs by interaction with the surrounding color 
field. This results finally in the third stage of a physical bound state. Both the second 
and third stages are non-perturbative in nature. 

On a fundamental theoretical level, color neutralization is not yet fully under- 
stood, but there are several models, color singlet [19], color octet [20], and color 
evaporation [5T]. The color evaporation model provides a particularly simple and 
experimentally well-supported phenomenological approach. In the evaporation pro- 
cess, the cc can either combine with light quarks to form open charm mesons {D 
and D) or bind with each other to form a hidden charm (charmonium) state. A 
fixed fraction of the subthreshold cc production is used in charmonium production. 
The basic quantity in this picture is the total sub-threshold charm cross section Sec, 
obtained by integrating the perturbative cc production cross section a over the mass 
window from 2mc to 2m_D. Since at high energy, the dominant part of Sec comes 
from gluon fusion (Fig. I1.10|l , we can write 



with gp{x) and gt{x) denoting the gluon densities and xi and X2 the fractional 
momenta of the gluons from projectile and target, respectively; a is the gg — > cc 
cross section. 

As mentioned, the basic assumption of the color evaporation model is that the 
production cross section for any particular charmonium state is a fixed fraction of 
the subthreshold charm cross section. 



PDF 




PDF 




(1.16) 
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a^is) ^ .h Scs{s) (1.17) 

where fi is an energy-independent constant to be determined empirically. It follows 
that the energy dependence of the production cross section for any charmonium state 
is predicted to be that of the perturbatively calculated sub-threshold charm cross 
section. As a further consequence the production ratios of different charmonium 
states 

^4 = Jl= constant (1.18) 
c^As) fj 

must be energy independent. Both these predictions have been compared in detail 
to charmonium hadro-production data over a wide range of energies (25]. They are 
found to be well supported, both in the energy dependence of the cross sections and 
in the constancy of the relative species abundances. 

Before turning to the topic of quarkonium production in hadron-nucleus colli- 
sions, let us consider the relevant time scales for the J/tp formation. 

The formation of a cc pair requires a time Tc- — l/2mc = 0.05 fm. The produced 
cc pair is in a color-octet state. To form a physical resonance state, it has to neutralize 
its color. The color-octet model [20] proposes that the color-octet cc combines with 
a soft coUinear gluon to from a color-singlet (cc — g) state. After a short relaxation 
time rg this pre-resonance (cc — g) turns into physical resonance by absorbing the 
accompanying gluon, with similar formation processes for the other resonances, such 
as Xc and ip' production. The color-octet model encounters difficulties if the coUinear 
gluons are treated perturbatively, indicating once more that color neutralization 
seems to require non-perturbative elements. However it does provide a conceptual 
basis for the evolution of the formation process. 



0.05 fm 0.25 fm 

hard | pre-resonance | resonance 

Tcc=l/2mc X8=l//2m|,Aqcd 



Fig. 1.11. Evolution of J/?/) formation in a hadron-hadron collision 



The color neutralization time Tg, of the pre-resonant state can be estimated by 
the lowest momentum possible for the confined gluons rs ~ {2mcAQCD)~^^^ — 0.25 
fm. The resulting scales of J/tp formation are illustrated in Fie: (|l.ll|l . The formation 
time for the actual physical ground state J/tp is presumably somewhat larger than 
Tg, although rj/^ — Ts , the heavy c quarks move non-relativistically. For the larger 
higher excited states , the formation time will then be correspondingly still larger. 

1.4.2 Quarkonium production in pA and AA collisions 

Let us now turn to nuclear collisions. Both in p-A and A-A interactions there will 
be pre-resonant absorption in nuclear matter. In nucleus-nucleus collisions, however. 
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there can be in addition a substantial amount of a produced "secondary medium", 
and testing this medium is in fact our main objective. 

The creation of the medium and production of the probe lead to two distinct 
formation scales. In p-A collision there is no formation time for the medium, so 
that such collisions provide a tool to probe charmonium production, evolution and 
absorption in confined matter. 

Nuclear effects can arise in all the evolution stages of J production, and a 
number of different phenomena have to be taken into account. 

• The presence of other nucleons in the nucleus can modify the initial state parton 
distribution functions, which enter in the perturbative cc production process, as 
shown in Fig. 11.101 This can lead to a decrease (shadowing) or to an increase 
(antishadowing) of the production rate. 

• Once it is produced, the cc pair in its evolution will traverse the nuclear matter; 
it can suffer absorption both in the pre-resonance as well as in the resonance 
stage, caused by successive interactions with the target nucleons. 



Hence J/t/j production in a nuclear medium is modified as compared to hadronic 
collisions. The modification occurs before any QGP formation and is thus indepen- 
dent of the effects due to a deconfined medium having free quarks and gluons. If we 
want to use J/ij) production and its suppression in a nuclear collision as a poten- 
tial signature of the QGP formation, all normal nuclear effects must first be taken 
into account. Only then can charmonium suppression serve as a probe to test the 
confining status of the produced "secondary medium" in nuclear collisions. 

So an essential question is how to account for the effects of the nuclear medium 
initially present on the production. The basis for this, both in theory and in ex- 
periment, is the measurement of dilepton, open charm and charmonium production 
in p-A or d-A collisions. These collisions thus provide a crucial tool to understand 
quarkonium production in nuclear collisions. 

The procedure to be used for such studies is the following: 

• The initial state parton distribution functions in nuclear matter are determined 
by open charm and dilepton production in p-A/d-A collisions in the relevant 
kinematic region. 



A 



A 




Fig. 1.12. J I'll) production in a nuclear medium. 
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• The Glauber model then is used to determine the pre-resonance absorption of 
the J/t}) and ip' by the target nucleon in p-A/d-A in the relevant kinematic 
region. 

It is thus clear that p-A or d-A coUision experiments are an absolutely essential tool 
for the analyis of quarkonium production in nuclear collisions. 

1.4.3 Sequential quEirkonium suppression 

There is a further important and, as it turns out, crucial feature observed in 
J/i/' hadroproduction. The J/ip actually measured in hadron-hadron collisions are 
not all directly produced 15 charmonium states; rather, they have three distinct ori- 
gins. About 60% of them are indeed directly produced IS charmonium states, but 
the rest are feed-down from higher excited states. About 30% come from the decay 
Xc(l-P) — > J/V' + anything and the remaining 10% from i/)'(2S') — > J/tp -\- anything. 
In both cases, the decay widths of the involved higher excited states are extremely 
small (less than 1 MeV), so that their lifetimes are very long and the decay occurs 
long after the interaction. The presence of any medium in nuclear collisions would 
therefore affect these excited states themselves and not their products, and we had 
seen above that excited states are dissociated before the ground state is. This has a 
direct consequence on the nature of J/'j/j suppression by deconfinement. In a ther- 
mal QCD medium, we should expect that with increasing temperature or energy 
density, first the J/ij) originating from %p' decay and then those from Xc decay will 
disappear. Only a considerably higher temperature would be able to remove the di- 
rectly produced J/ijjs,. Such a stepwise onset of suppression with specified threshold 
temperatures is perhaps the most characteristic feature predicted for charmonium 
production in nuclear collisions. It is illustrated schematically in Fig. 11.131 where 
we have defined the J/ip survival probability to be unity if the production rate suf- 
fers only the estimated normal nuclear suppression. The generic suppression pattern 
shown here will of course be softened by nuclear profile effects, impact parameter 
uncertainties, etc. On the other hand, this could be partially compensated if there 
is a discontinuous onset of deconfinement as a function of energy density of the 
medium. 



(2S)(1P) (IS) 



e(2S)8(lP) E(1S) 
Energy Density 




Fig. 1.13. Sequential J/^l> suppression 
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We had seen above how to calculate the quarkonium dissociation points which 
specify the temperature and thus also the energy density of the medium, thereby 
serving as a QGP thermometer. Potential model studies based on the heavy quark 
internal energy, as well as direct lattice QCD calculations gave as dissociation tem- 
peratures T ~ 1.1 Tc for tj}' and Xc and T > 1.5 — 2 Tc for J/tp. If this is correct, 
then the direct J/?/>(lS) survives up to about e > 10 — 20 GeVfm"^. Consequently, all 
anomalous suppression observed at SPS and RHIC must be due to the dissociation 
of higher excited states Xc and tp' [SS]- The suppression onset for this is predicted 
to lie around e ~ 1 GeV/fm'^, and once these are gone, only the unaffected J/ip 
(IS) production remains. Hence the J/tp survival probability (once normal nuclear 
effects are taken into account) should be same for central Au-Au at RHIC as for 
central Pb-Pb collisions at SPS. 

1.4.4 Charmonium regeneration 

In this section we want to investigate the possibility that the medium produced 
in high energy nuclear collision differs from the deconfined state of matter studied 
in finite temperature QCD. The basic idea here is that nuclear collisions initially 
produce more than the thermally expected charm and this excess, if it survives, 
may lead to a new form of combinatorial charmonium production at hadronization 
[Ml [251 [26]. 

In the QGP argumentation, a crucial aspect was that the charmonia, once disso- 
ciated, can not be recreated at the hadronization stage, because of the extremely low 
thermal abundance of charm quarks in an equilibrium QGP. The thermal production 
rate for a cc pair relative to a pair of light quarks is 

— ~ exp -(2mc - 2mq/Tc ) (1.19) 
~ exp(-2mc/rc) ~ 3.5 x 10"'^, 

with rric = 1.3 GeV for charm quark mass and Tc = 175 MeV for the transition tem- 
perature. The initial charm production in high energy hadronic collisions is, however, 
a hard non-thermal process, and the resulting rates calculated from perturbative 
QCD are considerably larger than the thermal rate. Moreover, in AA interactions 
the resulting c/c production rate grows with the number of binary collision A'coii, 
while the light quark production rate, being soft process, grow as the number of 
participants nucleons, i.e much slower. At high collision energies, the initial charm 
abundance in AA collisions is thus very much higher than the thermal value. Now 
the question is, what happens to this in course of the collision evolution? 

The basic assumption of the regeneration approach is that the initial charm 
excess is maintained throughout subsequent evolution, i.e, the initial chemical non- 
equilibrium will persist up to the hadronization point. In charmonium hadropro- 
duction, J/ip are formed because some of the cc pairs produced in a given collision 
form the corresponding bound state. In a collective medium formed through the su- 
perposition of many nucleon-nucleon (NN) collisions, such as a quark-gluon plasma, 
a c quark from one NN collision can in principle also bind with a c from another 
NN collision ("new" pairs) to create a J/tp. This pairing provides a "exogamous" 
charmonium production mechanism, in which the c and c in a charmonium state 
have different parents, in contrast to "endogamous" production in pp collision. At 
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sufficiently high energies this can lead to an enhancement in J/t/j production in AA 
collisions compared to a scaled pp rates [241 1251 126] . provided the overall charm 
density is sufficiently high at hadronization and provided the binding probability 
between charm quarks from different sources is large enough. 

Whether or not such enhancement becomes significant depends on two factors. 
On one hand, the initial charm oversaturation must be preserved so that the to- 
tal charm abundance is non-thermal. On the other hand it is necessary that the 
recombination between charm quarks from different parents to charmonium {J/ip) 
is strong enough. Here it is generally assumed that the final hadronization occurs 
according to the available phase space. Thus the number of statistically recombined 
J/ip has the form Nj/^ ~ N^^, growing quadratically in the number of cc pairs. 
This implies that the hidden to open charm ratio, e.g., Nj/^/Nd ~ Ncc/Nh in- 
creases with energy, in contrast to the energy independent form obtained for the 
fully equilibrated QGP, or to the decrease predicted by color screening. The pre- 
diction for J/ip production by regeneration is compared in Fig. 11.4.41 to that from 
sequential suppression. 



Fig. 1.14. Statistical J/i/) regeneration vs. sequential J/tj} suppression 



1.5 Conclusion 

Statistical QCD predicts the existence of a new state of nuclear matter, the quark- 
gluon plasma (QGP), at very high temperatures and/or densities. This medium, in 
contrast to hadronic matter, is capable of dissociating quarkonia, so that quarko- 
nium suppression may be taken as a sign of QGP formation in nuclear collisions [3]. 
Furthermore, different quarkonia dissociate at different temperatures; the dissocia- 
tion pattern thus serves as a "thermometer" for the QGP. It is therefore important 
to obtain precise predictions for these dissociation points, and for this, one can turn 
to either of two approaches, potential models or lattice studies. The former have the 
problem that the results are dependent on the type of potential chosen, while the 
latter so far suffer from the fact that lattice spacing and statistics limits the resolu- 
tion of peak widths in the spectrum. It is also not easy to identify the continuum 
region of the spectrum on the lattice. 




Energy Density 
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But what happens in case of relativistic nuclear collisions in the laboratory? If 
there is no regeneration of the dissociated charmonia, J/i/; remain as an external 
probe, and the sequential suppression pattern of the J/i/' can then serve as a tool 
to determine the energy density and the temperature of the produced medium. 
On the other hand, if there is J/^j) production through statistical combination of c 
and c from different collisions, leading to an overall J/ip enhancement, this would 
clearly indicate the thermalization of the produced medium on a pre-hadronic level. 
However, charmonia could then no longer serve as a thermometer to charaterize the 
primordial medium. Data from LHC, soon to come, will certainly play a decisive 
role in settling the issue. 
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